On the stability of U(l) spin liquids in two dimensions 



O 

o 

(N 

> 

O 

O 



c3 



a 

o 

> ■ 

in : 

r- . 
^t- ■ 
o ■ 
^t- : 
o . 
■ 

c3 ' 



a 

o 

o 



X 
S3 



Michael Hermele, 1 T. Senthil, 2 Matthew P. A. Fisher, 3 Patrick A. Lee, 2 Naoto Nagaosa, 4 and Xiao-Gang Wen 2 

1 Department of Physics, University of California, Santa Barbara, California 93106 
Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
3 Kavli Institute for Theoretical Physics, University of California, Santa Barbara, California 93106 
J * CREST, Department of Applied Physics, University of Tokyo, Tokyo 113, Japan 

(Dated: February 2, 2008) 

We establish that spin liquids described in terms of gapless fermionic (Dirac) spinons and gapless 
17(1) gauge fluctuations can be stable in two dimensions, at least when the physical SU(2) spin 
symmetry is generalized to SU(N). Equivalently, we show that compact QED3 has a deconfined 
phase for a large number of fermion fields, in the sense that monopole fluctuations can be irrelevant 
at low energies. A precise characterization is provided by an emergent global topological U(l) 
symmetry corresponding to the conservation of gauge flux. Beginning with an SU(N) generalization 
of the 5 = 1/2 square lattice Heisenberg antiferromagnet, we consider the 7r-flux spin liquid and, 
via a systematic analysis of all operators, show that there are no relevant perturbations (in the 
renormalization group sense) about the large- iV spin liquid fixed point, which is thus a stable phase. 
We provide a further illustration of this conclusion with an approximate renormalization group 
calculation that treats the gapless fermions and the monopoles on an equal footing. This approach 
directly points out some of the flaws in the erroneous "screening" argument for the relevance of 
monopoles in compact QED3. 



I. INTRODUCTION 

Phases of strongly correlated electrons in two dimen- 
sions exhibit a remarkable array of unusual and inter- 
esting behavior. Quantum spin liquids 1 - (Mott insula- 
tors with no broken symmetries) are a particularly exotic 
class of states that may play a critical role in our under- 
standing of several interesting materials. While much 
theoretical progress has been made in understanding the 
universal properties of such states, many basic questions 
remain unanswered. One of these concerns the stabil- 
ity of spin liquids described at low energies by gapless 
fermionic S = 1/2 spinons and a fluctuating £7(1) gauge 
field. We focus here on the cases where the spinons have a 
linear dispersion in the vicinity of discrete Fermi points 
and can be described as Dirac fermions. One example 
of such a phase, the staggered flux (sF) spin liquid pro- 
posed in Ref . 0, may play a key role in understanding the 
underdoped cuprate superconductors 3,4 ' 5,6 . 

In this class of U(l) spin liquids the fermions and 
gauge fluctuations interact strongly at low energies, and 
there are no well-defined low-energy quasiparticles£i&i24. 
Therefore the issue of whether these states can be stable 
is in general difficult to resolve, but becomes tractable, 
for example, when the physical SU (2) spin is generalized 
to SU(N) and N is taken large^. The problem can be for- 
mulated in terms of N species of Dirac fermions coupled 
to a compact {7(1) gauge field in two spatial dimensions, 
often referred to as compact QED 3 . The gauge fluctu- 
ations are suppressed at large N and calculation within 
the framework of a l/N expansion is possible, but even in 
this case there has been significant controversy. Based on 
an analysis of the monopole fluctuations^ and the sym- 
metries in the sF stated, it was argued that the sF spin 
liquid is a stable phase and the gapless spin excitations 
are protected^ at least in the large N limit. However, 



Refs. ITITIlI argue that, due to a screening effect in 
3D, the monopole fluctuations always result in the con- 
finement of the C/(l) gauge field, thus destabilizing the 
sF state and its close relatives. In this paper we resolve 
the ongoing debate and argue rather rigorously that the 
compact QED3 problem is deconfined for large N and 
that at least some {7(1) spin liquids are indeed stable in 
two dimensions. Whether similar deconfinement also ob- 
tains in physically important SU(2) spin models is left 
open for future work. 

These spin liquids generally possess at low energies 
much higher symmetry than the microscopic spin mod- 
els in which they arise. Of particular importance is an 
emergent topological global {7(1) that is associated in the 
gauge theory description with conservation of the gauge 
flux. In a compact {7(1) gauge theory flux is only con- 
served modulo 27r due to the presence of monopole fluc- 
tuations, so this emergent symmetry is equivalent to the 
irrelevance of monopoles at low energy. When monopoles 
are relevant it is believed that confinement is inevitable, 
so this symmetry provides precise meaning to the notion 
of deconfinement of spinons in the spin liquid. We em- 
phasize that our use of the term 'deconfinement' does 
not imply that the spinons are to be thought of as free at 
low energies. Roughly speaking, they are instead as free 
as possible given their strong interaction with the gauge 
fluctuations. A close parallel exists with the deconfined 
quantum critical points 14,15 that have been studied re- 
cently where an emergent topological global {7(1) also 
characterizes the low energy fixed point. Indeed the spin 
liquids discussed in the present paper may be viewed as 
deconfined quantum critical phases. 

We shall consider an SU(N) generalization of the 
5 = 1/2 Heisenberg model on the square lattice2ii£*i£. 
We define the model in terms of slave fermions f Ta 
(a = 1, . . . , N) transforming in the fundamental rep- 
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resentation of SU(N). Choosing TV even and impos- 
ing the local constraint /jj/ a = N/2 puts the spin in 
the antisymmetric self-conjugate representation. Defin- 
ing S af) = flfp - jjfiaf3f\f~f, the Hamiltonian is 

<rr') 

where the sum is over pairs of nearest-neighbor sites and 
the exchange interaction is assumed positive, J > 0. For 
N = 2 this reduces to the familiar 5 = 1/2 Heisenberg 
antiferromagnet, with global SU(2) spin rotation sym- 
metry. 

This model can be solved in the N — > oo limit, where 
slave fermion mean-field theory becomes exact. The cor- 
responding mean field Hamiltonian is quadratic and de- 
scribes the hopping of fermionic spinons with a hopping 
matrix element that is determined self-consistently. It 
is known that among mean field solutions that preserve 
all lattice symmetries the energetically favored one has 
a flux of 7r through every plaquette of the square lat- 
tice that is seen by the spinons. This state, known as 
the 7r-flux (ttF) state, may also be made the global min- 
imum energy mean field state upon addition of suitable 
biquadratic^ or ring exchanged terms. Here we are only 
concerned with stability of the 7rF state, which is a uni- 
versal feature independent of such microscopic param- 
eters. At the mean field level the spinon dispersion has 
two distinct gapless Fermi points. Linearizing the disper- 
sion in the vicinity of these points provides a description 
of the low energy physics in terms of a continuum theory 
of 2N species of 2-component Dirac fermions. 

The crucial question is the fate of this picture upon in- 
cluding fluctuations beyond the mean field. In the large- 
N spin model the primary effect of fluctuations is to in- 
duce a coupling of the fermionic spinons to a compact 
(7(1) gauge field. In (2 + 1) dimensions the compactness 
means that there are pointlike instantons (also known as 
monopoles) in space-time. At any such monopole event 
the total gauge flux associated with the (7(1) gauge field 
changes by an integer multiple of 2n. In pure (7(1) gauge 
theories (i.e. without any dynamical spinon fields) such 
instantons always proliferate and lead to confinement 19 . 
This then raises the question of whether the spinons of 
the mean field ttF state can escape confinement once the 
coupling to the compact (7(1) gauge field is included. In 
contrast to some recent studies jiLiSiiii we will argue that 
in the large N limit spinon deconfinement will indeed 
obtain. 

Consider first a description of the low-energy physics 
of the 7rF state that ignores the compactness of the (7(1) 
gauge field. The appropriate effective field theory con- 
tains 2N flavors of two-component Dirac fermions tpj 
(j=l,2,..., 2iV) minimally coupled to a noncompact (7(1) 
gauge field (with ^ = x,y,r). After rescaling spatial 
coordinates to set the Dirac fermion velocity to unity, the 



imaginary time action S — J drd 2 xC is given simply by 

C = tpjj^id^ + ia^j + — ^(e^A^aA) 2 , (2) 

one 

where the 7 M are 2x2 matrices satisfying {7^,7"} = 
25^. We have included an explicit Maxwell term for 
the gauge field, which will be generated by integrating 
out short wavelength modes. Earlier work has argued 
that this theory is critical and flows to a scale-invariant 
fixed point for N sufficiently large. The question then is 
whether properly accounting for the compactness of the 
gauge field, and including all other perturbations allowed 
by the microscopic symmetries, destabilizes this scale- 
invariant fixed point - in other words are monopoles (or 
some other allowed operator) a relevant perturbation in 
the renormalization group sense? 

The issue of stability to monopole fluctuations is sub- 
tle, and before attacking the problem with fermions it will 
be useful to briefly recap Polyakov's argument for con- 
finement in pure compact (7(1) gauge theory. By analogy 
with ordinary electrostatics, it is clear that monopoles 
will have a 1/r interaction in space-time coming from 
the Maxwell term in the above Lagrangian. If we let in- 
teger Qi represent the monopole number sitting on the 
sites of a 3D (say cubic) lattice, the effective action for 
the monopole gas will be given by, 

s m = \Y. QiQjVfc - pj ) + Qi > ( 3 ) 

with V(r) — e 2 n /r (for large r), where we have defined 
a "magnetic charge," e rn — l/2e. In the second term 
above, s c represents the action cost for the monopole 
cores, which will depend on short distance physics. Since 
a single monopole costs a finite action, due to entropic 
effects one expects them to always proliferate and be in 
a "plasma" phase. This can be simply established by de- 
coupling the monopole interaction term using a Hubbard- 
Stratanovich field Xi, an d then tracing out the monopoles 
on each site. For E c >> 1 one thereby obtains, 

S sg = \ J ^V- l {q)\x{q)\* - zYtKxixi), (4) 

where V(q) — Ane^/ q 2 and x(q) are Fourier transforms 
of V(r) and Xi, respectively, and z = 2e~ Sc . Back in real 
space the sine-Gordon action can be written in terms 
of a Lagrangian density, S sg — J d 3 rC sg , which in the 
continuum limit is simply 

C *9 = ^~( 5 ^) 2 - zcos (x), (5) 

where we have defined a coupling constant g = e 2 n . In a 
Hamiltonian picture, the operator e lx adds 2tt magnetic 
flux, and can be thought of as a magnetic-flux creation 
operator. In space-time it creates a monopole, so that 
the coupling z is correctly interpreted as the monopole 
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fugacity. Implementing a simple momentum-space renor- 
malization group (RG) perturbatively up to second order 
in z gives the RG flow equations, 



dz , gA 
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dt 
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with £ the usual logarithmic length rescaling and A a 
high momentum ultraviolet cutoff. Here c is a cutoff- 
dependent positive dimensionless constant. The com- 
bination A m = gA/n is the effective scaling dimen- 
sion of the monopole creation operator e* x , and rele- 
vance/irrelevance of z depends as usual on whether the 
scaling dimension is larger/smaller than the space-time 
dimenson D — 3. Since the coupling g scales to zero, 
the monopole fugacity clearly grows, indicating a pro- 
liferation of monopoles on long length scales. The field 
X gets pinned in the minimum of the cosine potential 
(a "smooth" phase), which dominates over the gradient 
term and generates a mass for \. This corresponds to a 
spatially short-ranged effective interaction potential, so 
the monopole gas is clearly in a plasma phase that can 
screen effectively. Due to the proliferation of "magneti- 
cally" charged monopoles, this corresponds to the phase 
of the U(l) gauge theory which confines the "electric 
charges" - in this case the fermionic spinons. 

We now turn to the effects of the gapless fermions on 
the issue of monopole proliferation. Existing arguments 
for and against stability of the monopole- free critical the- 
ory begin with a kind of random-phase approximation 
(RPA), where one integrates out the fermions to obtain 
an effective action for the gauge field. Denoting by a* 
the transverse part of the gauge field, this action has the 
following highly schematic form: 

ScffM = N (J AlqlM 2 + J (d 3 q) 3 -^\a t \ 

q (8) 
As N — > oo one has at ~ l/VN for fluctuations with 
finite action, and the nonlinear terms are apparently 
higher-order in 1/N than the leading Gaussian term; 
the action is therefore truncated at Gaussian order. For 
this quadratic theory, the effective interaction between 
monopoles is very long-ranged, Vrpa(7) ~ ln(7"), growing 
logarithmically in the space-time separation. Moreover, 
the action of a single monopole configuration diverges 
logarithmically in the infrared (the system size, say). The 
coefficient of the logarithm is proportional to N and very 
naively can be equated with the scaling dimension A m 
of the monopole creation operator (however, see later). 
The argument for the stability of the irF phase is then 
based on two observations, (a) Since the (apparent) scal- 
ing dimension of the monopole creation operator in the 
large N limit is (much) greater than the space-time di- 
mension, monopole fluctuations are irrelevant^iSS, and a 



+ 0{a & t ) 



naive perturbation theory that includes the monopoles 
should be free of infrared divergences, (b) The global 
lattice and spin rotation symmetries of the underlying 
spin model impose constraints on the continuum (non- 
compact) theory of the 7rF state, precluding, for exam- 
ple, a mass term for the fermions. This is a projective 
symmetry described by a projective symmetry groupAS. 
If this symmetry is unbroken, it protects the gaplessness 
of the spinons and the U(l) gauge bosons against all per- 
turbative fluctuations within the monopole-free sector of 
the theoryi2i2i. 

More recent work has, however, argued that this rea- 
soning is too naive and raised important and troubling 
questions. Indeed, beginning with the same quadratic 
gauge field action obtained within RPA by integrating 
out the fermions, it has been argued that the "bare" 
interaction between monopoles, Vrpa(?") ~ ln(r), will 
be screened by fluctuating monopole- antimonopole pairs 
present at finite monopole density. Specifically, screening 
was argued to reduce the logarithmic interaction to a 1/r 
form - the same potential that arises from the Maxwell 
term alone in the absence of any fermionsiiA^. Then 
Polyakov's argument (shown above) was invoked to con- 
clude that monopoles would always proliferate and lead 
to spinon confinement. 

This argument was explicitly demonstrated in the re- 
cent RG treatment of Herbut and Seradjehi^, who ef- 
fectively treated a gas of monopoles interacting via a 
pairwise logarithmic interaction, as in Eq. @ except 
with V(r) — > Vrpa(^). They made explicit use of the 
duality transformation to arrive at a sine-Gordon the- 
ory as above, but with an anomalous kernel. They ar- 
gued that the presence of the monopole fugacity term 
will lead to a self-energy correction at second order, 
V^\q) - (<?)+£(<?), with a self-energy £(<?) ~ z 2 q 2 
at small q. The effective "screened" interaction between 
monopole pairs is then simply 
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l + S( 9 )^RP A (g)' 
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which takes an intuitive "RPA- like" form with S(g) the 
monopole density-density correlation function. Since 
Vkpa{q) ~ q 3 , the S(g)T4i,p A (g) term dominates in 
the denominator, giving V c s{q) ~ I/S(g) cx 1/q 2 or 
V e g(r) ~ 1/r. Following Polyakov's original argument, 
the 1/r interaction is further screened to become short 
ranged, signaling a proliferation of monopole events and 
confinement. It appears that regardless of the value of A?", 
the proliferation of monopoles destabilizes the 7rF state 
and leads to spinon confinement. 

A first hint at the fallacy of this argument is that in 
the presence of monopoles the usual RPA treatment of 
the gauge interaction obtained by integrating out the 
fermions to quadratic order in a M is not correct, even 
in the large- N limftiK*. This is because while the gauge 
fluctuations are indeed suppressed at large- A", they are 
suppressed with respect to some classical background con- 
figuration. Monopoles are drastic perturbations to the 
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uniform — background, so in principle one needs to 
do a different fermion functional integral for every con- 
figuration of monopoles. Another way to say this is that, 
in the effective action Eq. JSJ), the finite- action N — > oo 
fluctuations about a static background of monopoles are 
not controlled only by the Gaussian term. This happens 
precisely because this action is written as an expansion 
about the wrong classical background. The same situ- 
ation occurs in the large- iV CP(iV) model at its critical 
poin1j2£ and related models^; these are similar to the 
model of interest in this paper, but have bosonic mat- 
ter fields. In all these cases, the monopoles are described 
by a strongly-coupled non-Gaussian theory with a rather 
specific structure of multi-monopole interactions. There 
is thus no reason to believe that a treatment based on 
a model as a generic Coulomb gas of monopoles with 
just pairwise logarithmic interactions is legitimate. Ar- 
guments and results valid for such a generic gas could po- 
tentially be modified for the rather special monopole gas 
that correctly obtains in the large-A limit in the present 
problem. 

But the fallacious monopole screening argument does 
raise an important point. To discuss the stability of 
Eq. J5J) to monopoles it is not enough to simply show 
that the scaling dimension of the monopole operator is 
larger than the space-time dimension. It is equally im- 
portant to ensure that as the monopole fugacity renor- 
malizes towards zero, the monopoles do not induce other 
operators in the monopole- free sector which are relevant. 
Indeed the screening occurring in the above sine-Gordon 
theory (i.e. the self-energy contribution) may be viewed 
as providing an explicit example where precisely this sort 
of thing happens. 

Thus the key to demonstrate stability of the irF state 
described by Eq. is to show that there are simply 
no relevant perturbations with or without the inclusion 
of monopoles. In this paper, we shall give a systematic 
analysis of all operators in the 7rF state and show that 
there are indeed no relevant perturbations, at least in the 
large N limit. 

An important closely related issue is that, in most ex- 
isting discussions of this problem, scant attention has 
been paid to the constraints imposed by the microscopic 
symmetries of the original lattice spin model on the ef- 
fective theory Eq. J5J). Apart from global SU(N) spin 
rotation these symmetries include translations, rotations, 
parity and time-reversal, as well as "charge conjugation," 
which takes the lattice fermion fields f a to f\. (For 
SU(2) spins charge conjugation is a subgroup of spin 
rotations, corresponding to a rotation by n around the 
y— axis in spin space, but for N > 2 is an independent 
discrete global symmetry of the spin Hamiltonian.) As 
emphasized in Refi ll ft the transformations of the fermion 
fields and the gauge field under these symmetries define 
the projective symmetry group associated with the slave 
fcrmion-gauge formulation and are crucial in prohibiting 
a class of potentially relevant operators in the continuum 
(monopole-free) theory. As mentioned above, without 



the projective symmetry, there is nothing to prevent the 
spinon mass term (a relevant perturbation) from being 
induced by any generic perturbation to the theory (be it 
inclusion of monopoles or other less drastic perturbations 
such as four- fermion interactions). 

II. GENERAL ARGUMENT 

We now present the details of our argument. The 
physics of the 7r-flux state and the associated fluctuations 
is encapsulated in the lattice gauge theory Hamiltonian 

HV{1) = \ E <Zr> (10) 
(rr') 

+ ••• 

where the ellipsis represents perturbations consistent 
with the symmetries. In the second term, the first 
sum is over sites of the A sublattice and the second is 
over nearest neighbors of r. Here a rr ' is a 27r-periodic 
vector potential living on the nearest-neighbor bonds, 
and e rr / is its canonically conjugate integer-valued elec- 
tric field. We must also specify the gauge constraint 

J2r'~> r e "" + fafa = lj tne nrst t erm is a lattice di- 
vergence of the electric field. With this choice Eq. I|10|) 
reduces exactly to the spin model in the limit h/t — ► oo 
with J ~ t 2 /h. The apparent breaking of lattice sym- 
metry is a gauge artifact, which is taken care of by re- 
quiring the lattice symmetry transformations to act on 
the spinons with additional gauge transformations, thus 
specifying the projective symmetry group of this spin 
liquid 1 — (see Appendix |A"|) . This information allows us 
to verify that the gauge theory Hamiltonian Eq. (|10|l has 
the exact same global symmetry group as the original 
spin Hamiltonian, and to determine what additional per- 
turbations are allowed. 

The continuum theory is obtained by first setting 
= and solving for the band structure of Eq. i|l(J|) 
(this reproduces the N = oo mean field state). Choosing 
a four-site unit cell, the spinon dispersion has nodes at 
(7r/2,7r/2) in the reduced Brillouin zone k x ,k y £ [0, n). 
For each component of SU (N) spin, the linear dispersion 
about the nodal point is described by two 2-component 
massless Dirac fermions. Then with N flavors of lattice 
spinons, one has a total of 2N flavors of 2-component 
Dirac fermions. Gauge fluctuations can be added to 
this continuum theory (as in Eq. J3J)) so long as we 
recognize that the gauge field is compact and allow for 
appropriate monopole configurations. It is straightfor- 
ward to determine the action of the microscopic sym- 
metries on the continuum fields; this is outlined in Ap- 
pendix^] This shows that all possible quadratic mass 
terms (i.e. fermion bilinears with no derivatives) are for- 
bidden. Moreover, the velocities associated with the dis- 
persion around the nodes are all required to be the same. 
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Thus the quadratic part of the Dirac action describes 2N 
flavors of 2-component Dirac fermions with full SU (2N) 
symmetry. Terms that break the SU(2N) symmetry are 
of quartic or higher order in the fermions, and as often 
happens the quadratic part has higher symmetry than 
the full action. 

It is important to observe that if monopoles are ig- 
nored, the theory in Eq. J5J) has in addition an extra 
topological global U(l) symmetry, as discussed in SecQ] 
This corresponds to the conservation of the gauge flux 
through any surface spanning the system in space-time. 
A charge-Q monopole event changes the flux by 2irQ and 
spoils this conservation law. 

Existing results in the literature show that the 
highly symmetric continuum theory Eq. (J2J flows to 
a conformally-invariant critical fixed point at large-iV. 
Here we examine all perturbations to this fixed point al- 
lowed in the present problem and argue that they are 
irrelevant. It is convenient to group the perturbations 
into two classes: operators that do not change the flux, 
and those that do. The former contains all perturbations 
that are allowed in the absence of monopoles. The lat- 
ter contains the monopoles, and their composites with 
polynomials in the fermion fields. 

Results from the 1/N expansion strongly suggest that 
the SU (2iV)-symmetric noncompact theory Eq. J3J flows 
to a critical fixed point for N sufficiently large£*&. The 
infinite- N theory is manifestly scale- invariant, with a 
photon propagator proportional to l/\q\ at small momen- 
tum. The effective expansion parameter for all diagrams 
in the 1/N expansion is then dimensionless and one ex- 
pects that only logarithmic divergences will occur; these 
contribute to nontrivial anomalous dimensions in princi- 
ple calculable order-by-order in 1/N. It should also be 
noted that a fermion mass cannot be generated pertur- 
batively in 1/N. (See Appendix [B] for a more detailed 
discussion on the validity of the 1/N expansion at large 
but finite N.) The criticality of the theory is further 
supported by a simple RG analysis in 3 — e dimensions^, 
which finds a nontrivial fixed point at 0(e) that is pre- 
sumably smoothly connected to the large-iV fixed point 
in d = 2. Furthermore, recent large-scale numerical stud- 
ies find no evidence for any symmetry breaking down to 
the smallest value of N simulated (N = 2)24. This con- 
clusion agrees with many approximate analyses of this 
model, finding SU(2N) symmetry breaking only below 
some critical jy<i25i2£i2£, contradicting early claims that 
the symmetry is always brokenSLSLSI. 

We next consider adding flux-preserving perturbations 
to Eq. that break its global symmetries down to 
those of the original lattice model. As emphasized earlier 
the microscopic symmetries forbid all possible quadratic 
fermion mass terms, as well as any velocity anisotropy. 
The leading such perturbations will therefore be four- 
fermion terms. At large-iV all 4-fermion operators have 
scaling dimension A — 0(1/N), and are hence expected to 
be irrelevant. If some such operator had instead been rel- 
evant, it could have led to spontaneous breakdown of the 



SU(2N) symmetry. Thus we expect that at least for N 
large the general non-compact theory (i.e. with only the 
global symmetries of the lattice model Eq. I|10|)) flows to 
an SU(2N )-symmetric conformally invariant fixed point. 

Finally we consider the monopole operators. Follow- 
ing Ref. |^ it will be extremely convenient to adopt a 
powerful point of view familiar to conformal field the- 
ory aficionados (for a pedagogical review accessible to 
condensed matter theorists see Ref. l32|) : In any D- 
dimensional conformally invariant theory, there is a one- 
to-one mapping between local operators and quantum 
states of the same theory quantized on the surface of 
the unit sphere S D . Furthermore, energy eigenstates 
on the sphere correspond to eigenoperators of RG scale 
transformations. The scaling dimension of such an oper- 
ator is equal to the energy of the corresponding quantum 
state. Ref. {3l| used this point of view to calculate very 
simply the scaling dimension and other quantum num- 
bers of monopole operators for the theory in Eq. J2J at 
leading order in 1/N. 

Consider a general local operator that changes the 
flux by an amount 2ttQ. This corresponds to a state 
on the sphere in a sector with a total magnetic flux of 
2ttQ through the surface. The energy of the state (and 
hence the scaling dimension of the corresponding oper- 
ator) is bounded below by the ground state energy in 
the same sector. In the N — > oo limit the gauge fluctu- 
ations are completely suppressed and the magnetic flux 
can be treated as a static, uniform background. Fur- 
thermore, the fermions do not interact in this limit, so 
the problem is reduced to finding the ground state en- 
ergy of 2N free Dirac fermions on the sphere in a back- 
ground magnetic field. As shown explicitly in Ref. 
(and is physically reasonable) the ground state energy 
in each such sector with non-zero flux is O(N). Con- 
sequently all flux changing operators have a scaling di- 
mension bounded below by a number of O(N) and are 
irrelevant. Moreover, as shown in Ref. 0, the scaling di- 
mension of a charge Q monopole is not equal to Q 2 times 
the scaling dimension of the charge Q = 1 monopole, as 
is the case within the Gaussian RPA treatment in terms 
of a gas of monopoles with a simple pairwise logarith- 
mic interaction. This makes clear that the conclusion in 
the monopole "screening" argument that the logarithmic 
interaction is screened down to a 1/r form rests on the 
flawed assumption of pairwise monopole interactions. 

We conclude that the effective theory in Eq. [21 flows to 
a conformally invariant fixed point with no relevant per- 
turbations - with or without the inclusion of monopoles. 
Consequently, the 7r-flux U(l) spin liquid state survives 
as a stable deconfined gapless critical phase, at least at 
large- N. 



III. RG ANALYSIS 

These considerations can be illustrated by the follow- 
ing (approximate) RG calculation which further helps 
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clarify the flaws in the RPA screening argument in 
Refs.]0 and^J We implement an approximate RG cal- 
culation following the approach used in Ref.0. However, 
our calculation pays attention to two important points: 
(a) we assume that the fermion mass counter term is not 
allowed (since it is forbidden by the projective symme- 
try in the irF state) ; (b) we let the infrared cutoff length 
scales for the fermions (Lf) and the monopoles (L m ) ap- 
proach infinity with a fixed ratio, say Lf/L m = 1. This 
is reasonable since both the monopole and the fermion 
sectors involve long-wavelength and gapless degrees of 
freedom. The choice in Ref. 0, with Lf taken to in- 
finity before L m , is not justified because it is based on 
an erroneous adiabatic approximation which requires the 
fermions to be rapidly varying variables compared with 
the monopoles. 

In the absence of gapless fermions, the RG flow equa- 
tions in the monopole sector in terms of the "running" 
monopole fugacity, z(£), and the running "magnetic 
charge," g(t) = e^(£), with I = ln(L m ), are already given 
in Eqs. JJjJ and J7J. Next consider the monopole-free sec- 
tor with gapless fermions coupled to the non-compact 
gauge field with "electric charge" e (with Lagrangian as 
given in Eq.(3J). An approximate RG perturbative in 
e can readily be implemented on this theory, which is 
simply noncompact QED3. This is accomplished by sim- 
ply ignoring all nonlinear terms in the action (except of 
course the minimal coupling vertex); in an exact treat- 
ment these would have nonzero coefficients at the fixed 
point of interest. Up to order e 6 the RG flow equations 
take the simple form, 



8^] 
di 



(4-D)e 2 -4tt?7— e 4 , 



(11) 



with t = ln(L^) and A an ultraviolet high momentum 
cutoff. Here rj > is a dimcnsionless number. The 
first term is present because the electric charge is not 
dimensionless in three space-time dimensions (D = 3), 
while the second term comes from a one-loop fermion 
bubble screening the Coulomb interaction and reducing 
the effective electric charge, as is familiar from four- 
dimensional QED. For large N this flow equation has 
a perturbatively accessible stable fixed point (with cou- 
pling (e*) 2 = A/ArjN) which describes the critical 7r— flux 
state. 

At this level of approximation, where we keep only 
the quadratic Maxwell term in the gauge action, the in- 
verse relation e m = l/2e is retained under the RG flows. 
Therefore we can deduce the effects of the fermions on 
the monopoles by replacing e 2 — > l/4g in the electric 
charge flow Eq. Ijlll) above, which gives, 



dg 



■K7]N 

IT' 



(12) 



The first term represents the engineering dimension of 
the magnetic charge in 3D, while the second comes from 
the screening of the fermions, which reduces the electric 
charge and thereby increases the magnetic charge. 



To arrive at a full set of RG flow equations incorpo- 
rating both the monopoles and the fermions, we must 
add to the right hand side of Eq. lfT2"|l the contribution to 
the magnetic charge coming from screening by monopole- 
antimonopole pairs, given explicitly in Eq. (|7|). The full 
set of coupled RG flow equations for the monopole fugac- 
ity and (magnetic) charge then take the form, 



dz . 5A 
at 7r 



dg 
dl 



-fJ 



cg 3 z 2 



Ttr/N 



(13) 



(14) 



If z(l) scales to zero in the I — > 00 limit, the mag- 
netic charge approaches a fixed point value, 5(00) = 
g* = irr]N/A. In this case the scaling dimension of the 
monopole creation operator becomes, A TO = g*A/n = 
r/N. For large N, A m is much greater than D = 3 so 
that z indeed does scale to zero consistent with the orig- 
inal assumption. Thus, provided 



Nrj > 3, 



(15) 



this approximate RG scheme predicts that the monopole 
fugacity scales to zero and the (magnetic) charge ap- 
proaches a fixed point value. The U(l) gauge field is 
not confining at such a fixed point, which corresponds to 
the 7rF phase with gapless spin excitations. 

One can see that the result Eq. I|15|) is the same as 
that obtained by Ioffe-Larkin 9 without considering the 
monopole screening effect. Their analysis corresponds to 
the single relation in Eq. I|13|l only, assuming g = g* ~ N. 
This is eventually justified in our coupled Eqs. i|T3|) 
and dJ. The arguments in Refs. ITllH2l focus on the — z 2 
term on the right side of Eq. I|14l) which represents screen- 
ing of the magnetic charge by monopolc-antimonopole 
pairs, a term which scales to zero in the present treat- 
ment. 

We emphasize again that the RG calculation above is 
approximate (even in the large- N limit). In particular 
the coupling of the gauge field to the fermions is treated 
in perturbation theory and all non-linear terms in the 
action are ignored. Nevertheless, it illustrates the failure 
of the monopole screening argument in this context. 



IV. DISCUSSION 

The low energy fixed point that describes the irF spin 
liquid has a global SU (2 A) symmetry and an extra global 
topological U(l) symmetry. The latter is a precise con- 
sequence of the asymptotic irrelevance of monopoles and 
corresponds to conservation of gauge flux. It thus pro- 
vides precise meaning to the notion of deconfmcment of 
the U(l) gauge field. These results may be expected to 
generalize to other states that also have Dirac fermions 
coupled to U(l) gauge fields in two dimensions such as 
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the staggered flux spin liquid. Our analysis, however, is 
controlled only at large- N, which simply provides a limit 
where these issues can be reliably addressed; future work 
will be required to determine whether similar deconfinc- 
ment obtains in models with real SU(2) spins. 

We remark that in three dimensions the stability prob- 
lem is essentially trivial even for SU(2) spin. In this case 
the large- N limit does not play a role; for example, it is 
known that U(l) spin liquids can exist as stable phases 
in d = 3 even if the spinons are completely gapped. It 
is very likely possible to have a stable three-dimensional 
liquid of SU{2) spins described by massless Dirac spinons 
interacting with an emergent U(l) gauge field - one need 
only show that all possible fermion mass terms are for- 
bidden by the microscopic symmetries. In this case the 
effective theory is QED in four space-time dimensions, 
and at low energies the spinons and photons interact only 
weakly. Because this theory is under good control and 
detailed predictions can be made, it may be very fruit- 
ful to search for possible realizations of such a phase in 
experiments and numerical simulations. 

We also note that issues similar to that resolved in 
this paper arise in the problem of a Fermi surface of 
spinons coupled to a compact U(X) gauge field in two 
dimensionsi^^. Since this problem has an even higher 
density of low lying excitations compared with Dirac 
fcrmions, the notion of integrating out all the fermions 
and considering monopoles in the resulting action trun- 
cated to quadratic order is even more suspect. We expect 
that a correct treatment will likely show that monopoles 
are irrelevant for large N in this case as well. Of course, 
here the monopole-free theory is likely vulnerable to 
Fermi surface instabilities. It is hence unlikely to sur- 
vive as the ground state but could perhaps be as stable 
as an ordinary Fermi liquid. 
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APPENDIX A: SYMMETRIES AND 
CONTINUUM FIELDS 

For completeness, we provide here a discussion of the 
continuum limit of the 7rF mean field state and the action 
of the microscopic symmetries on the continuum Dirac 
fields. The starting point is the mean field irF Hamilto- 
nian 

«* = -*££ [(*+(-i) (r *- ri) )/rU'«+k. C.]. 

(Al) 

Note that this is simply Eq. (|10f> with gauge fluctuations 
completely suppressed (a rr ' =0). 

It is convenient to work with a four-site unit cell la- 
beled by (R, i), with R = 2n x x + 2n y y and r(R, i) — 
R + Vj , where 



(A2) 



The spinon operator at the site (R, i) is denoted fjn a . It 
is a trivial exercise to go to momentum space and solve 
Eq. IjAlfl : in the reduced Brillouin zone k x ,k y £ [0, n) ap- 
propriate for this unit cell one finds gapless Fermi points 
at Qo = (7r/2,7r/2). Near this point the dispersion can 
be described by 27V 2-component Dirac fermions. It is 
convenient to denote these by tp^a C^) ■ Here a = 1 , 2 
and a = 1, . . . , N are the SU(2N) flavor indices (a is 
simply the SU(N) spin index). Also, A = 1,2 labels the 
two components of each spinor (this is often suppressed) . 
These fields are related to the lattice spinons as follows: 




^e^- R (/ RlQ + / R3Q ) (A3) 
lQ °- R (/R 2a - /R4«) (A4) 
e iQo - R (/R 2a + fn4a) (A5) 



2V2£ 

_ e ~iir/4 



2V2£ 



-iir/4 



2-J21 



-e lQoR (/Ri a -/R 3 a), (A6) 



where £ is the lattice spacing. 

In momentum space the continuum Hamiltonian takes 
the form 



d 2 q 



aa V My 



(A7) 



where we have chosen units to set the velocity to unity, 
and t 1 are the usual Pauli matrices acting in the 2- 
component Dirac "spin" space. Here we use the following 
rotated coordinates: 



32 



1 
1 



{-q x + %)■ 



(A8) 
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We now simply quote the action of the microscopic 
symmetries on the lattice spinons and the resulting trans- 
formations for the continuum fields. In order to keep Ti.^ 
invariant, in some cases the spinons transform with an 
additional U(l) gauge transformation; this is the hall- 
mark of projective symmetry and the projective sym- 
metry group, and has important consequences for the 
action of the symmetries on the continuum fields. We 
often include an extra uniform gauge transformation 
fra — > e 8 * fra to (slightly) simplify the form of the con- 
tinuum transformation laws. It is convenient to adopt a 
four-component notation, defining: 



*a = 



^2a 



(A9) 



We can represent matrices acting in this four-component 
space by the tensor products of Pauli matrices t 1 ^ . The 
t % matrices act in the Dirac spin space, while the /j. 1 act 
in the flavor space connecting the upper and lower 2- 
component spinors of \&. In the S'[/(2A r )-symmetric con- 
tinuum theory, the fx 1 generate the SU(2) subgroup of 
SU(2N) consisting only of emergent symmetries; these 
are in some sense continuous extensions of the discrete 
lattice symmetries. 

x- translations. Translations by one lattice site in the 
x-direction act on the spinons as follows: 



fr, 



ex P(-yCr)/r+x,a, 



(A10) 



where £ r = 0, 1, 2, 3 when the coordinates of r are (even, 
even), (odd, even), (odd, odd) and (even, odd), respec- 
tively. The continuum transformation law is 



(All) 



Rotations. We choose to make a ir/2 counterclockwise 
rotation about the center of a plaquette; that is, we rotate 
about the point (x+y)/2, which gives r = (r x , r y ) — ► r' = 
(—r y + 1, r x ). The action on the spinon operators is: 



where 



fra * ^-rfr'ai 



■IreA 
-1 reB 



(A12) 



(A13) 



In the continuum we have: 

^(R)^exp(f(^))exp(^3)^ (R0 

Reflections. We consider a reflection of the form 
r = (r x ,r y ) — > r' = (— r x ,r y ), under which the spinons 
transform trivially 



fra * fr'a 

In the continuum this leads to: 



(A15) 



* a (R) - (* M 2 ) exp (I^!))* Q (R/) (A16) 



Charge conjugation. As discussed in Sec. [H] charge 
conjugation is distinct from the SU(N) spin rotation 
symmetry only for N > 2. The action on the spinons 
is 



f ra -> exp(yCr)/ 1 



t 

rev 



In the continuum: 
4rt 



(JT 1 )^ 1 )^ 



(A17) 

(A18) 
(A19) 



Time Reversal. Time reversal is an antiunitary opera- 
tion acting on the lattice spinons as follows: 



fra * f r/m 

In the continuum this becomes: 



(A20) 

(A21) 
(A22) 



It is at this point a simple exercise to show that 
these symmetries forbid all possible mass terms M^- = 
^t//V' 1 E , Q ,, as well as any velocity anisotropy. 



APPENDIX B: VALIDITY OF THE 1/N 
EXPANSION AT LARGE BUT FINITE N 



While our results are based on an analysis within the 
framework of the 1/N expansion, they should hold for 
large but finite N. In order to see this it is important to 
consider the precise connection between the 1 /N theory 
and actual finite- TV models. 

Consider the field theory Eq. (J2J for some finite value 
of N. We can formally integrate out the fermions to ar- 
rive at an effective action for the gauge field as in Eq. © 
- one can imagine using an appropriate lattice regular- 
ization in order to ensure that this step is well-defined. 
With the fermion fields gone, there is no obstacle to treat- 
ing iV as a continuous variable and doing perturbation 
theory in 1/N. The resulting family of field theories is 
manifestly non-local because we have obtained it by inte- 
grating out the gapless fermions; however, for integer N 
it is equivalent to the local field theory we started with. 
For non-integer N there is no reason to believe the the- 
ory is equivalent to any local theory. The key point is 
that near 1 /N = the points that do correspond to a lo- 
cal theory become arbitrarily closely spaced, so it should 
be a very good approximation to view the theories in a 
small, continuous interval near 1/N = as local. 

The above considerations mean that we can think of 
our expansion about the 1 /N — fixed point in the stan- 
dard framework of the Wilsonian renormalization group, 
which applies only to local theories. In this case a finite 
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value of 1/N is an exactly marginal perturbation, and 
we can calculate corrections to the scaling dimension of 
any operator as 1/N is tuned away from zero. The 1/N 
expansion itself gives only an asymptotic series for each 
of these scaling dimensions. However, if the asymptotic 
N — > oo result for a given scaling dimension A is greater 



than the space-time dimension D, it is straightforward 
to use the mathematical definition of asymptotic conver- 
gence to show that A > D for some sufficiently large but 
finite value of N, above which the corresponding operator 
is irrelevant. 
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